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ABSTRACT 

Recent observations have discovered a population of extended Lya sources, dubbed Lyman- 
alpha blobs (LABs), at high redshift z ~ 3 — 6.6. These LABs typically have a luminosity of 
L ~ io 42 ~ 44 erg s" 1 , and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 
kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger 
model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts 
by means of a combination of hydrodynamics simulations with three dimensional radiative transfer 
calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass 
range of 3 — 7 x 10 12 M© in the redshift range of z ~ 3 — 7, and we use the newly improved ART 2 
code to perform the radiative transfer calculations which couple multi-wavelength continuum, ioniza- 
tion of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling 
induced by gravitational interaction produce strong Lya emission from these merging galaxies. The 
Lya emission appears to be extended due to the extended distribution of sources and gas. During 
the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10 3 M Q yr -1 , our 
model produces Lya blobs with luminosity of L ~ 10 42 ~ 44 erg s -1 , and size of D ~ 10 — 20 kpc at 
z > 6 and D ~ 20 — 50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. 
Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant 
ones may be produced by mergers more massive than those in our model, or a combination of mergers 
and cold accretion from filaments on a large scale. 

Subject headings: galaxies: high-rcdshift - galaxies: formation - galaxies: evolution - radiative transfer 
- radiative feedback - hydrodynamics - methods: numerical 



1. INTRODUCTION 

Star formation in galaxies can produce strong Lya 
line emission via recombina tion in Hli region in the 
early phase of the ev olution dPartridge fc Peebles! 119671 ; 
ICharlot fc Fal [1993ft . Recent narrow-band deep imag- 
ing surveys using large-aperture telescopes have detected 
a large population of Lya emitting galaxies, generally 
called Lyman-alpha emitters (LAEs) when the equiv- 
alent width (EW) of the Lya line i s above 20 A , a t 
nearly all reds hift up to z = 8.6 (e.g..lCqwie fe Hulll99£ 
Hu et al.f i 998: Mal hotra fe Rho ads 2004; live et allboOt 



Gronwall et all l2007t lOuchi et alf 1201(1 iLehnert et all 
20101: iKashikawa et all 120111 : iShibuva et al.1 120121 1 It 

is suggested that mos t LAEs are compact, young star- 
forming galaxies fe.g.. iWestra et al.ll2005l iGawiser et~al 



120071: lAcauaviva et all 120121 but see iKornei et al.ll20iC 

for a different view). 

Meanwhile, a number of extended and bright Lya 
sources, which are called Lyman-alpha blobs (LABs), 
have also been discovered in th e redshift range o f 
z ~ 3 - 6.6 (IKeel et all H99l ISteidel et all 1200 



20111: iReuland et al.l 120031: iMatsuda et al l 12004 I201T 



Dev et all 120051: iSaito et all I2006L l2008t iNilsson et al 



20061: iSmith fc Jarvisl 120071: iPrescott et al.l ~l200__ 
Ouchi et al .1120091: 1 Yang et al T I2009L 120101 ). These LABs 
have a luminosity in the range of ~ 10 42-44 erg s -1 , 
a size from tens of kiloparsecs up to ~ 100 kpc, much 
more extended than the LAEs at the same redshift, and 

|yuhl9@psu.edu| 



some LABs do not seem to have any optical or infrared 
counterpart. 

The origin of the LABs remains a hot debate. One 
proposal is that they may form fro m gravitationa l cool- 
ing radiation from accreting gas (|Haiman et al.l 120001 : 
iFardal et al.l l200l"|) . motivated by the observations that 
some LABs lac k visible power source in both opt i- 
cal and infrared ([Nilsson et al.ll200l ISmith et al.|[200l . 
Recent hydrodynamics simulations show that galaxy 
evolution is accompanied by accretion of cold gas 
streams of T ~ 10 4 ~ 5 K, which penetrate deep inside 
dark matter haloes (|Katz et ail 120031: iKeres et "ail 1 20051 
20091: iBirnboim fc Dekell l2003t iDekel fe Birnboirri 12009: 



Ocvirk et al.ll2008tlBrooks et all 2009; D ekel et al.l200£ 



Such inflow of cold gas may produce strong and extensive 
Lya emission via collisional excitation, which may result 
in LABs if there is no dust absorpt i on (iDiikstra fc Loebl 
20091 iFaucher-Giguere etall [20091 iGoerdt et all 120101: 



Yaiima et al.ll2012cD . 



However, other observations shows that some LABs 
are associated wit h reg ular Lyman break galax- 
ies (jMatsuda et al.l 12004ft . sub-millimeter and in- 
frared sources which imply active star formation 

of SFR ~ 10 2 - 3 M<T) /yr (IChapman et al.1 120011 

iGeach et al.l 120051 120071). or active galactic nuclei 



(AGN ) and quasars dBunker et aLl)2003t iWeidinger et all 
2001 IBasu-Zvch fc Scharfl 120041: IGeach et all \mm 



Smith et all 120091 : iColbert et all 1201 lh . More recently 



observations showed as ymmetric, filamentar y structure 
in so me giant LABs (IMatsuda et all 120111: lErb et all 
120111: iRauch et all 120111 12012ft . which may result from 



cold accretion or interaction-triggered inflow. 

To explain the extended distribution of Lya emis- 
sion of LABs, stellar wind from supernovae feed- 
back has been suggested as an effective mec h anism 
(|Taniguchi fc Shioval l2000( ) . iMori fe UmelmiTal (p006l ) 
showed, by means of ultra-high resolution hydrodynam- 
ics simulations, that multiple supernovae feedbacks pro- 
duced extended bubble structure and the compressed gas 
shell created numerous Lya photons over extended re- 
gion. However, it is not clear that triggers the starburst 
or AGN activity. In addition, most of the theoretical 
works assume that all Lya photons can escape from the 
galaxy, but in reality, Lya photons can experience nu- 
merous scattering and be absorbed by dust. Moreover, 
galactic wind is highly s uppressed by gr avitation poten- 
tial at M ha i > 10 12 M Q (|Springelll200a ). which may not 
be effective in producing LABs. 

In this work, we propose a formation model of LABs 
from major mergers of gas-rich galaxies. In such 
a merging process, strong gravitational torques cre- 
ate highly condensed gas clumps in the nuclei region 
and tidal t ails, and trigger intense, global star forma- 
tion (e.g.. iHernouistl 119891: iBarnes fc Hernauistl 119921: 



Sanders fc Mirabell Il996t ISpringell 120051: IHopkins eTaLl 
20061) . At higher redshift (e.g., z > 2), the major merger 



rate is higher in de nse regions and the progenitors ar e 
more gas rich (e.g., iLi et all 120071 IHopkins et all 120101 ) . 
which may produce strong, extended Lya emission from 
both recombination of ioni zing photons and co llisional 
excitation of hydrogen gas (lYaiima et al.H2012allc1 lbh. 

To test this model, we investigate the Lya property of 
interacting galaxies by combining hydrodynamical simu- 
lations with three dimensional radiative transfer (RT) 
calculations. The simulations follow the evolution of 
galaxy mergers of different mass and redshift, and the 
RT calculations uses the newly improved code ART 2 by 
lYaiima et all (|2012aD . The ART 2 code couples multi- 
wavelength continuum, Lya line, and ionization of hy- 
drogen, which is critical to study the Lya and multi-band 
properties of galaxies. 

The paper is organized as follows. We describe the 
galaxy simulations in §2, and the method of radiative 
transfer calculations in §3. In §4, we present the results 
of Lya properties, which include the photon escape frac- 
tion, emergent Lya luminosity, the EW, the line profile, 
Lya surface brightness and the size. We discuss in §5 
the effect of wind and AGN on the Lya properties, and 
summarize in §6. 



2. GALAXY MODEL 

In order to determine the role mergers play in pro- 
ducing LABs, we perform a set of hydrodynamics simu- 
lations of binary, equal-mass, major mergers of galax- 
ies .^_Tliegalaxyis constructed based on the model 
of I Mo et all (|1998l ), which consists of a dark matter 
halo, a disk of gas and stars, and a s eed black hole of 
10 5 M ( n, using a well-t e sted method (IHernauistl [19931: 
Springel fc White! Il999l : ISpringell 120001: ISpringel et alJ 
2005bJ). We follow the procedures of ILi et all (|2007l ) to 
generate the galaxy progenitors, the properties of which, 



including the virial mass M v j r , virial radius i? v ir and 
halo concentration C v i r , arc scaled appropriately with 
redshift. 
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where G is the gravitational constant, and Mq ~ 8 x 
1O 12 /i _1 M0 is the linear collapse mass at the present 
epoch. 

The density pr ofile of the dark matter halo follows a 
Hernquist profile (|Hernq uist 1993), scaled to match tha t 
found in cosmo logical simulati ons (|Navarro et al.|[l~997T) . 
as described in ISpringel (|2005l ). 
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where a is a parameter that relates the IHernauistl (|1990f ) 
profile parameters to the appropriate NFW halo scale 
length R s and concentration C v - U - (C v - lr = R V1I /R S ) 7 



a = iW2[ln(l + C v ir) - C vir /(1 + C vir )]. 



(6) 



The exponenti al disk of st ars a nd gas are then con- 
structe d as in IHernquis"tl (|1993D and ISpringel et al.l 
(|2005bl) . We assume a baryon fraction of /b = 0.17 for 
these high-redshift galaxies based on the seventh year 
Wilkinson Microwave Anisotropy Probe data (WMAP7, 
iKomatsu et aLll2011h . The gas fraction of each progen- 
itor is extrapolated from t he results of semi-ana lytical 
models of galaxy formation (jSomerville et al.ll2001l ) , with 
100% gas disks at z > 10, 90% at 10 > z > 6, and 60% 
at 6 > z > 4. 

The simulations include dark matter, gasdynam- 
ics, star formation, black hole growth, and feedback 
processes, and are performed using the parallel, N- 
body/Smoothed Particle Hydrodynamics (SPH) code 
GADGET -3, which is an impr o ved version of that de- 
scribed in ISpringel etaD (l200ih : ISpringell (f2005h . GAD- 
GET implements the entropy-con serving formulation of 
SPH (TSpringcl fc Hernauistl 2002T) with adapt ive parti- 
cle smoothing, as in iHernquist fc Katzl (|1989l ). Radia- 
tive cooling and heating processes are calculated assum- 
ing collisional io nization equilibrium (|Katz et al.l 119961 : 
iDave et"aT1ll999ft . Star formation is modeled in a multi- 
phase interstellar medium (IS M), with a rate t hat follows 
the S chmidt-Kennicutt Law (jSchmidtl [T959; Kcn nicuttl 
119981 ). Feedback from supernovae is captured through a 
multi-phase model of the IS M by an effective equation o f 
state for star-forming gas (ISpringel fc HernquistJ I2003T) . 
The UV background model of lHaardt fc Madaiif (|1996f ) 
is used. A self-regulated black hole model is used, in 
which a sphe r ical Bondi-Hoy l e-Lytt l eton gas accretion 
(IBondil H95l IBondi fc HovH [19441 IHovle fc Lvttletor] 
Il941| ) with an upper limit of Eddington rate is assumed, 
and the feedback is in the f orm of thermal e nergy i njected 
to the gas, as described in ISpringel et al.l (|2005al) . 

We also include a model of ga lactic wind driven by 
stellar feedback, as introduced by ISpringel fc Hernquist! 
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TABLE 1 

Simulation runs with different parameters and feedback models. 



Runs 


A/tot (M Q ) a 


Zimt 


JVdm c 


N h d 


BH 


Wind 


/gas 6 


z9A 


3 x 10 12 


9 


6 x 10 5 


4 x 10 5 


on 


on 


0.98 


z6A 


5 x 10 12 


6 


6 x 10 5 


4 x 10 5 


on 


on 


0.9 


z4.5A 


7 x 10 12 


4.5 


6 x 10 5 


4 x 10 5 


on 


on 


0.6 


z4.5B 


7 x 10 12 


4.5 


6 x 10 5 


4 x 10 5 


off 


on 


0.6 


z4.5C ! 


7 x 10 12 


4.5 


6 x 10 5 


4 x 10 5 


on 


off 


0.6 



a Total mass of the merging system. 

Initial redshift of the simulation. 
c Total number of dark matter particles 
^ Total number of star and gas particles 
c Gas mass fraction to dark matter normal 



(|2003f ). We adopt a constant wind velocity of t> w ind = 
484 km s _1 , a mass- loss rate that is twice of the star for- 
mation rate, and an energy efficiency of unity such that 
the wind carries 100% of the supernova energy. The wind 
direction is anisotropical, preferentially perpendicular to 
the galactic disk. This wind model causes an outflow 
of gas, transporting energy, matter and metals out of 
the galactic disk in proporti on to the star formation rate 
(jSpringel fc Hernquistl 12003 ) . 

Table 1 lists the simulations performed in this work 
and the related physical and numerical parameters. The 
three fiducial runs, z9A, z6A, and z4.5A are set up to 
produce galaxy mergers at redshift z ~ 6, 4, and 3, re- 
spectively, for direct comparison with observations. The 
progenitors of these mergers belong to the massive end 
at the given redshift. All these three simulations in- 
clude thermal feedback from supernovae and black holes, 
and galactic wind. In order to investigate the effects of 
feedback on the Lya properties of these merging sys- 
tems, we also run two additional simulations, z4.5B and 
z4.5C, without black hole or wind, respectively. We 
adopt the cosmological parameters from the WMAP7 re- 
sults: H = 72kms" 1 Mpc- 1 (h = 0.72), Q M = 0.26, 
Q A = 0.74, Q b = 0.044, a$ = 0.80, and n s = 0.96. 

3. RADIATIVE TRANSFER 

The RT calculations are performed using the 3D Monte 
Carlo RT code, All-wavelength Radiative Transfer with 
Adaptive Refinement Tree (ART 2 ), as recently developed 
by lYaiima et all (l2012al). A RT 2 was improved over the 
original version of iLi et all ()2008l ). and features three 
essential modules: continuum emission from X-ray to 
radio, Lya emission from both recombination and col- 
lisional excitation, and ionization of neutral hydrogen. 
The coupling of these three modules, together with an 
adaptive refinement grid, enables a self-consistent and 
accurate calculation of the Lya properties, which depend 
strongly on the UV continuum, ionization structure, and 
dust content of the object. Moreover, it efficiently pro- 
duces multi-wavelength properties, such as the spectral 
energy distribution and images, for direct comparison 
with multi-band observations. The detail ed implementa - 
tion s of the ART 2 code ar e described in ILi et all ()2008l ) 
and lYajima et all (|2012a| ) . Here we focus on the Lya 
calculations and briefly outline the process. 

The Lya emission is generated by two major mecha- 
nisms: recombination of ionizing photons and collisional 
excitation of hydrogen gas. In the recombination process, 
we consider ionization of neutral hydrogen by ionizing 



ized by WMAP-7 year data 

radiation from stars, active galactic nucleus (AGN), and 
UV background (UVB), as well as by collisions by high- 
tcmpcrature gas. The ionized hydrogen atoms then re- 
combine and create Lya photons via the state transition 
2P — > IS. The Lya cmissivity from the recombination is 



f a aBhv a n e nHu, 



(7) 



where ae is the case B recombination coefficient, and 
f a is the average number of Lya photons produced 
per case B recombin ation. Here we use as derived in 
IHui fc Gnedml (fl997l ). Since the temperature depen- 
dence of f n is not strong, = 0. 68 is assumed every- 
where (|Osterbrock fc Ferlandll2006l ) . The product hv a is 
the energy of a Lya photon, 10.2 eV. 

In the process of collisional excitation, high tempera- 
ture electrons can excite the quantum state of hydrogen 
gas by the collision. Due to the large Einstein A co- 
efficient, the hydrogen gas can occur de-excitation with 
the Lya emission. The Lya emissivity by the collisional 
excitation is estimated by 



oil r < 



Lya«o«HI, 



(8) 



is the collisional excitation coefficient, 

-17. 



3.7 x lQ- 17 exp(-h u n /kT)T- 1 / 2 ergs s" 1 cm 3 



where Ch ya 

Ci, Vfl = 3.7 x iu "exp(— /» /, 
(|Osterbrock fc Ferlandll200l . 

Once the ionization structure have been determined, 
we estimate the intrinsic Lya emissivity in each cell by 
the sum of above Lya emissivity, e a = e™ c + e£° u . 

In RT calculations, dust extinction from the ISM is in- 
cluded. The dust content is estimated according to the 
gas content and metallicity in each cell, which arc taken 
from the hydrodynamic simulation. The dust-to-gas ra- 
tio of the MW is used where the metallicity is of Solar 
abundance, and it is linearly interpolated for other metal- 
licity. We u se the stellar population synthesis model of 
GALAXEV (jBruzual fc Cha riot 2003]) to produce intrin- 
sic spectral energy distributions (SEDs) of stars for a grid 
of metallicity and age, and we us e a s imple, br o ken p ower 
law for the AGN (|Li et all 120081 ). A lSalpeterl (fl955h ini- 
tial mass function is used in our calculations. 

We apply ART 2 to the above hydrodynamic merger 
simulations. In our post-processing procedure, we first 
calculate the RT of ionizing photons (A < 912 A) and 
estimate the ionization fraction of the ISM. The result- 
ing ionization structure is then used to run the Lya RT 
to derive the emissivity. Our fiducial run is done with 
N p h = 10 5 photon packets for each ionizing and Lya 
compone nts, which was demon strated to show good con- 
vergence (|Yaiima et al .11201 2 al io!) . The highest refinement 
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of the adaptive grid corresponds to a cell size comparable 
to the spatial resolution of 20 pc in physical coordinate 
of the hydrodynamic simulation. 

4. Lya PROPERTIES 

The simulations z9A, z6A, and z4.5A produce mergers 
at redshift z ~ 6, 4 and 3, respectively. The high-rcdshift 
galaxy progenitors in our simulations are gas-rich and 
compact. The merging systems typically experience a se- 
ries of separation, close encounter, and final coalescence 
phases. During the merging event, strong gravitational 
interactions between the progenitors lead to tidal tails, 
strong shocks and efficient gas inflow that triggers large- 
scale starbursts and fuels rapid accretion onto the black 
holes. Meanwhile, feedback from both supernovae and 
black holes suppress both star formation and black hole 
accretion. The resulting Lya emission follows the evo- 
lution of the system, as illustrated in Figure [T] from the 
z4.5A run. 

At t = 0.1 Gyr after the simulation starts, the galaxy 
pair is undergoing the first close encounter. The Lya 
emission extends to a large distance and shows two strong 
peaks around the two interacting nuclei. As the pair coa- 
lesces at t <~ 0.48 — 0.5 Gyr, in contrast to a strong spike 
of starburst, the Lya emission dims and shrinks, and 
shows holes in the central region of the system, indicat- 
ing significant absorption. After that, the Lya emission 
continues to decline and its distribution becomes more 
compact in the merger remnant. 

Figure [2] shows the star formation history and the re- 
sulting Lya properties of the three fiducial simulations, 
z9A, z6A, and z4.5A, respectively. Higher redshift sys- 
tem merges faster due to smaller dynamical time scale. 
All simulations show intense star formation during the 
merging process, which in general declines with time ex- 
cept a strong peak at the time of the final coalescence, 
with rates ranging from a few to ~ 10 3 Mq yr . 

The majority of the Lya photons are absorbed by dust 
in these merging systems. The escape fractions of Lya 
photons, / e8c , are below 30%, and during the final coa- 
lescence phase, they fall to deep dips of / osc ~ 1%. This 
is caused by the sudden increase of absorption from the 
dust produced by the starburst, as well as a sudden in- 
crease of the scattering optical depth in the shocked and 
compressed gas regions. As a result, the Lya luminos- 
ity declines with time, in particular, the z6A run shows 
a steep drop at the coalescence, even its SFR reaches 
~3x 10 3 M Q yr -1 . The z4.5A run generally has higher 
luminosity, L-^ya ~ io 42 - 7-44 - 2 erg s _1 , the z9A having 
L-Lya ~ if) 41 - 7-43 - 8 erg s _1 , and z6A lies in between. 
The observed luminosity range of LABs in z ~ 3 — 6 
is ~ 10 43 — 10 44 erg s _1 , within the range of our simula- 
tions. 

The equivalent width (EW) of Lya line is defined as 
the ratio of Lya luminosity to the UV flux density at 
> 1216A, EW = Lj^ya/fw, where /uv is the mean flux 
from 1300 A to 1600 A. As shown in Figure H the EWs 
of z4.5A run shows ~ 25 - 220A, while z9A has ~ 40 - 
100A, and z6A - 50A. This range is in broad agreement 
with observation s of EWs of LABs, - 40 - 200A (e.g., 
iSaito et all 120011) . 

From the three-dimensional Lya map, we measure the 
size of the blob in a projected plane above a given flux 



limit. If the two galaxies are separated apart, we mea- 
sure the size of the bigger one only. Figure [3] gives the 
size evolution of the Lya blobs from the three fiducial 
simulations above 2.2 x 1 0~ 18 erg s" 1 cm~ 2 ar csec" 2 , 
the detection threshold of iMatsuda et all (|2004h . The 
z4.5A run has a size in the range of D ~ 15 — 50 kpc 
during the merging event, z6A in ~ 15 — 45 kpc, and 
z9A in ~ 10 — 20 kpc. These results are consistent with 
observations. For example, the most distant LAB de- 
tected at z = 6.6 bv lOuchi et al.l (j2009T ) has a luminosity 
Lhya = 3.9 x 10 43 erg s _1 and a size D ~ 17 kpc, in 
good agreement with that in the z9A simulation, which 
represents a merger event in z ~ 6. 

5. DISCUSSION 
5.1. Effects of Galactic Wind and AGN Feedback 

It has been suggested that stellar wind from super- 
novae or AGN feedback may produce the extended Lya 
distribution of LABs (e.g.. iTaniguchi k Shioval 120001: 
iMori k Umeliural 120061 ; lOuchi et all 120101 ) . To exam- 
ine the effects of feedback on the formation of LABs, we 
compare simulations with and without feedback. 

In our si mulations, we adopt the supe rnovae feedback 
model fro mlSpringel k Hernquistl (|2003l ). and AGN feed- 
back from lSpringel et al.l (|2005al T The feedback from su- 
pernovae includes both thermal and kinetic forms. The 
thermal energy is released to the surrounding ISM, which 
results in a multi-phase structure. The kinetic feedback 
is in form of an anisotropic galactic wind perpendicular 
to the galactic disk. The wind has a constant velocity 
of Wwind = 484 km s _1 , a mass- loss rate that is twice of 
the star formation rate, and an energy efficiency of unity 
such that the wind carries 100% of the supernova energy 
For the AGN feedback, about 5% of the radiated energy 
is injected to the gas thermally and isotropically, which 
would raise the gas temperature and cause an isotropic 
outflow. 

Figure [4] shows a comparison of the LAB luminos- 
ity and size among simulations z4.5A (with both super- 
novae and AGN feedback), z4.5B (no AGN), and z4.5C 
(no wind, but with thermal feedback from both super- 
novae and AGN). While the feedback does not seem to 
affect significantly the resulting Lya luminosity, the su- 
pernovae wind appears to make a remarkable difference 
in the Lya distribution. Without wind, the Lya is more 
concentrated, reducing the size by ~ 10%. However, the 
AGN feedback does not seem have a significant effect on 
the Lya distribution. 

5.2. Limitations of Our Model 

While our model can produce luminous (^Lya ~ 
10 43-44 erg s -i) &nd exten d c d (£) ~ 20 - 50 kpc) Lya 
sources at z ~ 3 by idealized, binary major mergers, they 
fall short to reprod uce the observed gian t LAB s with a 
size of - 100 kpc (jMatsuda et al.l 12004 120 111) . These 
giant LABs may form from mergers more massive than 
the ones we modeled, or from multiple mergers, or from 
a combination of mergers and accretion of cold gas along 
the filaments on a large scale. 

In a cold dark matter cosmology, galaxies form hier- 
archically, with small objects collapsing first and sub- 
sequently merging to form progressively more massive 
ones. In this picture, massive galaxies form in high over- 
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Fig. 1. — Evolution of the Lya map during the merging event, from the z4.5A simulation. The box size is 100 kpc in physical scale, and 
the color bar indicates the surface brightness of the Lya flux. 



densi ty regions where mergers occur frequently (|Li et al.l 
[2001 . The merging event may occur along the filamen- 
tary structure of intergalactic medium, and may be ac- 
companied by accretion of cold gas onto the galaxy from 
the filaments. In these cases, the Lya emission may have 
very extensive distribution due to scattering along the 
cold flows and filaments. 

To systematically investigate the formation of LABs, 
an ideal approach would be to combine multi- wavelength 
radiative transfer calculations, as those presented here, 
with full cosmological simulations which follow galaxy 
formation and evolution at different cosmic epochs. Such 
large-scale simulations, however, typically lack the high 
resolutions required to resolve the Lya emission. We will 
tackle this issue in future work. 

6. SUMMARY 

By combining hydrodynamical simulations of interact- 
ing galaxies and radiative transfer calculations, we inves- 
tigate the possibility of producing extended Lya blobs 
from binary, major mergers. We focus on three mergers 
taking place in the rcdshift range of z ~ 3 — 7 with a mass 
range of 3 - 7 x 10 12 M . We find that the highly shocked 
gas regions and intense star formation induced by the 
gravitational interaction produce strong Lya emission 
from these merging galaxies. The Lya emission appear to 
be extended due to the extended distribution of sources 



and gas. It has a luminosity of L ~ 10 erg s 1 
and a size of D ~ 10 - 20 kpc at z > 6, D ~ 20 - 50 
kpc at z ~ 3. These results are in good agreement with 
observations of typical Lya blobs from redshift z ~ 3 to 
6.6. Furthermore, we find that feedback from supernovae 
or AGN does not have a significant impact on the Lya 
flux, but the presence of galactic wind from supernovae 
appear to increase the size of the blobs by ~ 10%. 

We note that our model does not produce giant Lya 
blobs of ~ 100 kpc. This suggests that giant LABs may 
form from mergers more massive than the ones we pre- 
sented here, or from multiple mergers, or from a combi- 
nation of mergers and cold accretion along the filaments 
on a large scale. To fully account for the Lya sources 
in these scenarios, one needs to combine high- resolution 
cosmological hydrodynamics simulations of galaxy for- 
mation and evolution, with three-dimensional radiative 
transfer calculations, which we plan to investigate in a 
forthcoming paper. 
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Fig. 2. — Evolution of the physical and Lya properties during 
the galaxy mergers from the three fiducial simulations, z9A, z6A, 
and z4.5A, which produce an merging event at z ~ 6, 4, and 3, 
respectively. From top to bottom in clockwise direction are: the 
star formation rate (upper left panel), the emergent Lya luminosity 
(upper right), the equivalent width of Lya line (lower right), and 
the escape fraction of Lya photons (lower left). 
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